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Observations of the high energy sky, mainly with the INTEGRAL satelHte, have raised new ques- 
^ r"| tions about the formation and evolution of High Mass X-ray Binaries (HMXBs). The number of 

detected HMXBs of different types is now high enough to allow us to carry out a statistical analy- 
O sis of their distribution in the Milky Way. For the first time, we derive the distance and absorption 

^ of a sample of HMXBs using a Spectral Energy Distribution fitting procedure, and we examine 

^ the correlation with the distribution of Star Forming Complexes (SFCs) in the Galaxy. We show 

that HMXBs are clustered with SFCs with a typical cluster size of 0.3 ± 0.05 kpc and a charac- 
teristic distance between clusters of 1.7 ± 0.3 kpc. Furthermore, we present an investigation of 
the expected offset between the position of spiral arms and HMXBs, allowing us to constrain age 
and migration distance due to supernova kick for some sources. These new methods will allow 
CO us to assess the influence of the environment on these high energy objects with unprecedented 

CN reliability. 
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1. Introduction 

Thanks to the dedicated observations from RXTE and INTEGRAL, around 200 HMXBs are 
currently known in the Milky Way allowing us to focus on their distribution. Using RXTE data, 
[9] highlighted clear signatures of the spiral structure in the spatial distribution of HMXBs. In 
the same way, [7], [11] and [2] and [3] showed that HMXBs observed with INTEGRAL also seem 
to be associated with the spiral structure of the Galaxy. However, the HMXB positions, mostly 
derived from their X-ray luminosity, are not well constrained and highly uncertain in case of direct 
accretion as in HMXB. In order to overcome this caveat we present a novative approach allowing 
us to uniformly derive all HMXB positions. We study this distribution and the correlation with 
Star Forming Complexes (SFCs) observed in the Galaxy. Knowing the location of these sources, 
one can examine the composition of the environment at their birthplace. This study is necessary to 
better understand the formation and evolution of the whole population of HMXBs and primarily to 
state the role of the environment and binarity in the evolution of these high energy binary systems. 
For a complete overview of this study, the reader should consult [6]. 

2. Deriving the HMXB location within our Galaxy 

2.1 SED fitting procedure 

Using the Liu et al. catalogue ([10]), updated with literature, and the INTEGRAL source 
catalogue of [1], we retrieved the quiescent optical and NIR magnitudes (mostly from the 2MASS 
point source catalog), the spectral type and the luminosity class of each source from the literature. 
We then built a sample of 46 sources for which at least 4 optical/NIR magnitudes and the spectral 
type are known. For each source we build its optical/NIR SED and fit it with a blackbody model. 
Two parameters are fixed in the fitting procedure: the radius and the temperature of the companion 
star which dominates the optical and NIR flux. Two other parameters are left free: the extinction 
in V band Ay and the ratio R/D, whereas the extinction A^ is derived at each wavelength from [5] 
assuming Ry = 3.1. Knowing the radius R of the companion star, we then calculate the distance 
D in kpc. Results on distance and extinction determination compared with previously published 
results are given in [6]. We point out that median discrepancy in distance is ^ 17% and Ay is often 
very similar (median discrepancy of less than 7%). 

2.2 Uncertainties in the computed distance and extinction 

To compute accurate distance and extinction uncertainties, we estimated the magnitude un- 
certainties (retrieved from the literature), the uncertainties on the radius and temperature of the 
companion star, the degeneracy between several parameter values (based on the fitting procedure), 
the infrared excess of Be stars due to their circumstellar envelope and the error due to a different 
extinction law. For a complete analysis of these uncertainties see [6]. We present in Figure 1 the 
distribution of HMXBs in the Galaxy, obtained with our novative approach based on uniform dis- 
tance determination. The spiral arms model given by [12] is also presented. Figure 2 represents 
most of the studied HMXBs with the uncertainties on their location computed taking into account 
all the errors described above. The question then arises: is there a correlation between this distri- 
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bution of HMXBs and the distribution of Star Forming Complexes (SFCs) in the Milky Way (given 
by [12]), as it is expected from the short HMXB lifetime ? 
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Figure 1: Distribution of HMXBs (blues stars) and SFCs (green circles). The circle radius of SFCs repre- 
sents the excitation parameter value. The spiral arm model from [12] is also plotted and the red star at (8.5 ; 
0) represents the Sun position. 



3. Results: HMXB distribution and correlation with Star Forming Complexes 

To probe the correlation between HMXBs distribution and Star Forming Complexes, we use 
the method described in [6]. A correlation is detected with two characteristic scales: the typical 
cluster size of 0.3 kpc and the typical distance between clusters of 1.7 kpc with uncertainties of 
0.05 kpc and 0.3 kpc respectively. If we take into account the uncertainties in HMXB positions and 
in SFC positions (given in [12], median error of 0.25 kpc), the correlation still exists with the same 
cluster size and the same distance between clusters. [3] mention that HMXBs and OB complexes 
are clustered for a cluster size r < 1 kpc. This upper limit, obtained with a different method, is 
consistent with our results. 

4. Implication of the correlation on HMXB formation and evolution 

The distribution of HMXBs reflects the stellar formation that took place some tens of Myr ago, 
since they are not an instantaneous star formation rate (SFR) indicator as explained in [13]. Then, 
an offset between spiral arms (an indicator of the actual star formation) and HMXB distribution 
is expected. Indeed, since the spiral arm rotation velocity is different than the angular velocity 
of the stellar disk, we expect HMXB positions to be offset from the currently visible SFCs. This 
time lag was mentioned in [11] and [7] but a deeper investigation of this issue was not possible 
due to the limited sample of HMXBs at that time. [13] evaluated the offset for the galaxy M51. 
Here, we attempt to implement this formalism in the case of the Milky Way. One can assume the 
spiral arm density wave to rotate with the speed Q.p = 24 km s~^ kpc~^ (see e.g. [8]), in the same 
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direction as the stellar disk, which velocity curve fl(r) is assumed to be flat in the galactocentric 
distance range of interest, according to [4]. Then, to a first approximation, it is possible to derive 
the expected HMXB locations relative to the current position of the spiral arms in time T, i.e., the 
angular offset A0(r) given by A0(r) = {^(r) — Q.p)T. To estimate this displacement, we plot the 
expected positions of sources formed 10, 20, 40 Myr ago in the Galaxy map (Figure 2). For better 
visibility, here, we do not plot the expected positions of sources formed 60, 80 and 100 Myr ago. 
Even if Figure 2 does not inmiediately suggest by eye any substantial offset between the current 
spiral arm position and the expected position of HMXB s (which depends on the age of the sources) 
we would now like to quantitatively assess the presence of an offset. 

4.1 Existence of an offset between HMXBs and Galactic spiral arms 

To perform this study we must split the sample of HMXBs depending on the age of the sources. 
Two different samples are then created: one containing 4 supergiant stars and a second one con- 
taining 9 Be stars (the way these samples were created is explained in [6]). Then, we calculate 
the distance from each HMXB to the closest actual spiral arm (given in [12]). We follow exactly 
the same procedure to calculate the distance from each source to the closest expected position of 
sources formed 20, 40, 60, 80 and 100 Myr ago and we determine the mean value of the offsets 
(taking into account all the sources of the two samples). The sources are expected to be closest to 
one of the expected positions computed above than to the current spiral arms observed by [12]. We 
determine the l-<7 error bars to the distance between each source / and the closest point on the arm 
using the propagation of uncertainty formula. Taking into account all the sources of the sample, 
we cannot observe any significant variation of the offset with time before 60 Myr. This result was 
expected because the minimum peak of distance between each HMXB and the closest expected 
position should be clearly different whether we consider Be or supergiant stars. Also, by taking all 
spectral types into account, we tend to loose a part of the information except the mean age upper 
limit of 60 Myr highlighted by the plot. For both supergiant and Be samples, we observe a more 
significant increase of the offset after 60 Myr. This increase could be a signature of the expected 
offset between the HMXB positions and the current spiral arms. However, we must be cautious 
about this result especially because of the small number of HMXBs in the two samples. 

Moreover, the time interval during which HMXBs appear should translate the mass range of 
both stars of binary systems (see [7]). These results only enable us to state that this time interval is 
lower than 60 Myr on average for all stars. 

4.2 Deriving the age of HMXBs 

Following the method described above, one can compute a distance between each source and 
the different theoretical arms that correspond to the current predicted position of a source sample 
born 20, 40, 60, 80 or 100 Myr ago. In the previous section we studied a sample of Be and 
supergiant HMXBs, and this could also be applied to each source separately, to constrain the age 
and the potential migration of the system due to a supernova kick. 

It is possible to determine the rough age of these sources and a lower and upper limit of kick 
migration distances. If we consider only the sources studied separately here, we can derive a mean 
age of ^ 45 Myr for sgHMXBs and 51 Myr for BeHMXBs. This result seems consistent with the 
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Figure 2: Positions of HMXBs with error bars (zoom in the solar neighborhood). Color arms represent 
expected position of HMXBs with ages of 10 (blue), 20 (green) and 40 (red) Myr. 

distinct evolution timescales of these two kinds of systems but again, we must be cautious about 
this result because of the small sample of sources used in this calculation (4 supergiant and 9 Be 
systems). 

4.3 Constraints on migration from supernova kicks 

We focus on the sources selected before (4 sgHMXBs and 9 BeHMXBs). The distance be- 
tween the object and the closest expected position gives a kick value for each source and a mean 
value depending on the spectral type of the two samples. We derive a mean migration distance of 
0.1 1 kpc for BeHMXBs and of 0.10 kpc for sgHMXBs. However, we should be careful about these 
mean values because of large error bars and small samples. It is also important to underline that 
these derived values only represent a lower limit to the kick migration distance since we cannot 
take into account the migration distance on galactic latitude given by the kick. Then, we only get 
the projected migration distance on the Galactic plane. 

5. Conclusion 

Examining the distribution of HMXBs is of major interest in order to study in depth the for- 
mation and evolution of these high energy sources. However, HMXB locations are usually poorly 
constrained and largely dependent on the determination method. Here, we determine the location 
of a sample of HMXBs using an uniform approach: SED fitting of their distance and absorption. 
This method enables us to reveal a consistent picture of the HMXB distribution, following the spi- 
ral arm structure of the Galaxy. This study shows that HMXBs are clustered with SFCs with a 
cluster size of 0.3 ± 0.05 kpc and a distance between clusters of 1.7 zb 0.3 kpc. We go further 
by quantitatively assessing the offset between current spiral density wave position and expected 
HMXB positions. Moreover, for 4 sgHMXBs and 9 BeHMXBs, we are able to derive an age and a 
migration distance, giving constraints on the supernova explosion kick. 
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